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: Abstract 

[ We reconsider the DMO sum rule in light of our recent two-loop calcula- 

' tions of isospin and hypercharge vector and axialvector current propagators 

in chiral perturbation theory. A modified derivation valid to second order in 
the light quark masses is presented, and a phenomenological analysis yields 
! determination of a combination of finite counterterms occurring in the p 4 and 

p 6 chiral lagrangians. Suggestions are given for further study. 

> 

Working in the chiral limit of massless quarks, Weinberg was the first to derive spectral 
o3 ■ function sum rules involving the vector and axialvector currents. [|TJ There soon followed 
the announcement of two additional chiral sum rules, one for the pion electromagnetic mass 
difference and one (the DMO sum rule ||) involving an inverse moment of chiral spectral 
functions, both also derived in the chiral limit. A subsequent study of chiral sum rules in 
QCD determined that only the first Weinberg sum rule and the DMO sum rule survive the 
inclusion of quark mass effects. 0,0] About a decade later, the DMO sum rule was analyzed 
by applying the Gasser-Leutwyler analysis of chiral perturbation theory (ChPT) at one- 
loop order. |)p|| In this paper, we report on a re-analysis of the DMO sum rule based on 
our recent calculations of vector and axialvector current propagators at two-loop order in 
ChPT. H|0 1 



Derivation of DMO Sum Rule in Two-Loop ChPT 

Given the normalization of vector and axialvector octet flavour currents, 

J& = qjYq and J^ = q^f^q (< = !,... ,8) , (1) 
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one defines the isospin vector and axialvector current propagators respectively as 

A£ 3 V) =ijdf-x e** <0|T J£ 3 (0)) |0> (k = V, A) . (2) 

These have the tensorial decompositions 

^(^) = (^-^OnS(g 2 ) + ^g v ng(g 2 ) (k = V,A) , (3) 
and spectral functions obtained via imaginary parts, 

-ImU£tf) = W-fini®tf)+ffl$tf) (k = V,A) . (4) 

7T 

In the following we shall assume isospin symmetry, for which the spin-zero vector current 
contributions vanish (Ilyg = and py] = 0)- 

The asymptotic behavior (s — > oo) of the vector and axialvector spectral functions 
predicted by QCD is given as [§, @], JTT | 

(pS-pS-pS)W ~ °( s_1 ) • (5) 

The information in Eq. (|5[) can be used, together with analyticity and the corresponding 
asymptotic behavior of the polarization functions, to derive dispersion relations for the vector 
and axialvector polarization functions, 

(n«;> - ng - n<°») (9 >) = f ds ^ >' s » . (6) 

Upon evaluation at = 0, Eq. (|3|) yields the inverse- moment DMO (or 'Lio') sum 
rule 1, 

KS - nJ8 - nSS) (0) = f * (pl ' j " P "" P " >W • < 7 > 

J s 

For the purpose of phenomenological analysis, one wishes the spectral integral to contain 
only quantities which have already been measured. Therefore, we make use of the ChPT 
analysis of Ref. |12| to cast the DMO sum rule in the optimal form, 

n<'> - ft« + rf$g\ (0) + 4<( rt - 2L<?(„.) = fds , (8) 

where H[ ^ (/i) and (/i) are counterterms from the p 4 chiral lagrangian and nA3°' ) are finite 
polarization functions entering at two- loop order |T0|JT2]] . 

ChPT Contributions to DMO Sum Rule 

From Refs. [pl,[l0|1 , we obtain the following expressions for terms appearing on the left- 
hand-side (LHS) of Eq. (|), 

LHS = + G« = + G?> + Gf + + Ggj, , (9) 
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where and are respectively the ChPT one-loop and two-loop contributions and 

(2,) 

the final two (numerical) terms in G^u M arise from the so-called sunset amplitudes. All 
the two-loop quantities are new contributions which have not previously appeared in the 
literature. We shall discuss the counterterm contributions separately in a later section. 

Phenomenological Determination of DMO Spectral Integral 

One anticipates an accurate determination of the DMO spectral integral 

(pSJ-pS)W 



'DMO 



ds 



(11) 



insofar as the inverse factor of the squared-energy suppresses contributions at large s. 
The dominant contributions to Jdmo will arise from the 2-k and 3n resonances p(770) and 
ai(1260), with smaller contributions from n n > 4 multipion states, XX-multipion states, 
etc. One can access such information via studies of tau-lepton decay and e + e~ — > 2n, etc 



cross sections. 13 



An evaluation of Jdmo which employs as input the most recent r and e + e data available 
and which neglects contributions from s > m 2 yields ]14 



/dmo= (26.4 ±1.5) -10" 



(12) 



with an uncertainty of about 5.7%. In a modified analysis, one defines the laplace trans- 
formed quantity 



DMO 



(M 2 



s M 2 



ds e 



6M 6 24M 8 



(13) 



where Jdmo — liniM->oo Idmo{M 2 ) and the final two terms involve nonperturbative contri- 
butions. A fit to several M 2 -moments then yields fl~4| 



j DMO = (25.8 ± 0.3 ± 0.1) • 10~ 3 
with a claimed uncertainty of about 1.24%. 



(14) 
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Inputs to Counterterm Contributions 



Consider the contributions in Eq. (0) arising from counterterms. There are those from 



the chiral lagrangians of order p A (L^ and L^') and those from lagrangians of order 
(Qa S , R^ s , Qv S anc ^ -^v S )- All these counterterms are evaluated in MS renormalization.[] 
For notational consistency with the axial propagator counterterms, the vector propagator 
counterterms, defined earlier in Refs. |],[nj as Q and R, are denoted here respectively as 



(o)> 



Q>f and R ms 

The apparent dependence in Eq. (Q) on the scale \jl is illusory. The explicit scale depen- 
dence of the logarithmic terms is compensated by implicit scale dependence in the p A and p 6 
counterterms such that the LHS is independent of scale. In particular, the scale dependence 
of the p 6 counterterms is 
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as can be obtained from renormalization group equations ||. 

The collection of p A and p 6 counterterms originate from the renormalization procedure, 
and each must ultimately be determined from experiment. The p A counterterms and 



L^o are already known jjTjJ from one-loop analyses to accuracies respectively of about 10% 



and 13%. 



hf[M p ) = (6.9 ±0.7) • 10~ 3 
L$ (M p ) = -(5.5 ±0.7) • 10~ 3 



(from(r 2 )^) , 
(from 7r~ — > eVfPj) 



(16) 



We may employ these values for and in the two-loop quantities and G% of 
Eq. (^) because any error made is of still higher order. This is, however, not true for the 
L^o dependence in G^, for which consistency demands a two- loop evaluation. That is, 
suppose hypothetical^ that a physical amplitude existed which contained L^q , but no other 
counterterms to either one-loop or two-loop order. A phenomenological determination of 
lJ[o carried out in a two-loop analysis would yield a value somewhat different from that 
obtained in a one-loop analysis. We would expect this shift in L±q to be of the same order 
as the p 6 counterterms which comprise g£t in Eq. (0), and so we are not justified a priori 
in ignoring it. 

Of the p 6 counterterms, we have previously used spectral function sum rules to obtain 
the estimates, 



1 The reader is referred to Ref. [ 10 1 for a discussion of renormalization scheme dependence of the 
counterterms. 
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QT(M P ) 



(3.7 ±2.0) • 1(T 5 
(1.4 ±0.5) • 10~ 4 



(Ref. [15]) 
(Ref. 0) 



(17) 



where the result of Ref. \TU\ has been scaled down to fi = M p in the above. The relatively 



large uncertainties in each of these reflects the paucity of existing data in the hypercharge 
channel. It is also possible to obtain estimates of certain counterterms via the assumption 
of resonance exchange saturation JT7j. Thus, a reanalysis [|TI| of the paper by Jetter |19[] on 
77 — > 7T7T and i] — > ^77 yields 
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(18) 



Since the counterterm Qy s receives contributions in the resonance exchange approach only 
from scalar exchange, its magnitude (the sign is not fixed) is rendered small. The resonance 
couplings Cg and are defined in Appendix D.2 of Ref. 0], and the result is seen to agree 
nicely with our sum rule determination cited in Eq. (|i~7|). As for R^f s , although the contri- 
butions from low-lying resonances vanish, there will of course be continuum contributions. 
Even so, we disagree with the numerical estimates given in Ref. ||19|| , and the vanishing result 
for S I res implies that the interpretation of a large rj — > ^77 amplitude given in Ref. (19 



is not tenable. 



Analysis and Conclusions 

We are now ready to study the numerical implications of the DMO sum rule. To begin, let 
us temporarily avoid any use of counterterm inputs in order to get a result free of theoretical 
errors. Then using the conservative estimate of Eq. (0) for the DMO spectral integral Jdmo, 
we obtain from the DMO sum rule the numerical constraint 



Ml 

7T 

~F 2 ~ 



h 



Dvro _ \m p 



QT(M P )-QT(M P ) 

(2) 



M|±2M| r m 



F 2 



-< S (M p ) 



-0.0038 ±0.0004 



(19) 



where = G« + 4Lf ) and G {2) = G (2) - G^. The RHS of Eq. flU) will change in a 
known manner as the scale /z changes but the error bar, due exclusively to the uncertainty 
in the evaluation of -/dmo> Wl11 remain the same. The value of the RHS is seen to be rather 
smaller in magnitude than the one-loop value L^(M p ) ~ —0.0055. This is due in large part 
to sizeable chiral logarithms present in Gffi and G^. That is, the numerical contribution 
of Idmo to the RHS (equal to —0.0066) is reduced via cancellation with the one-loop term 

to -0.0054 and then via cancellation with the two-loop term G^ to the final value of 
—0.0038. The latter reduction is about a 30% effect. 

Can the contribution of L^{M p ) be disentangled from the mass corrections on the LHS 
of Eq. fllPp? The interest in this question is underlined by the large numerical effect, just 
mentioned, the chiral logarithms have on the LHS of the DMO sum rule. From Eq. (|i~9"D. 
one can already conclude that if L^(M p ) stays close to the one-loop determination cited 
in Eq. fll~6|) , the mass corrections involving the p 6 counterterms have to be both substantial 
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and negative at scale [i = M p . Conversely, if these mass corrections turn out to be small at 
jj, = M p , any two-loop analysis would need to yield a determination for L^(M p ) reduced by 
about 30% relative to its one-loop value. 

The information gathered in Eqs. (|T7| ) , (|T8l) is not by itself sufficient to conclusively answer 
this question. Nevertheless it is instructive to investigate in more detail the individual 
contributions to the LHS of the DMO sum rule. Using the estimates of Eq. (|17|) , the 
contribution of the Q-type counterterms is numerically found to be 

AM 2 r i 

- —f [QT(M P ) - QT(M P )} = (9.3 ± 5.0) • 10- 4 (20) 

where we have added errors on Q^ s and Qy S in quadrature. The Q-terms thus contribute 
approximately 3.5% to the value obtained phenomenologically for the RHS of the DMO 
sum rule. For i?y S we can employ the resonance saturation estimate of Eq. (|T8|). No 
corresponding estimate is available, however, for R^ s . A priori, such contributions proceed 
through exchange of scalar resonances as well, but the couplings of those resonances to two 
axialvector currents are not known experimentally. However, since both R-type counterterms 
are Zweig-rule suppressed, we expect these to be small compared to the Q-type counterterms. 

There is, however, a well-known drawback of such resonance saturation estimates of 
counterterms. The method does not specify the scale at which the counterterm has to be 
taken. The generally accepted procedure is to employ the resonance mass as the relevant 
scale and to use the variation of the counterterm in a range between, say, fi = 0.5 — > 1 GeV 
as an estimate of the uncertainty of the method. In the DMO sum rule, the counterterms 
turn out to vary strongly with scale. This can be seen either by using Eqs. (|15|) or be read 
off from Table 1 where we give the two- loop contributions Gffi , and Gnum f° r three 
different scales 

Table 1: Scale dependence 



li (GeV) 


0.50 


0.77 


1.00 




0.0018 


0.0023 


0.0026 




0.0001 


0.0042 


0.0065 


^NUM 


-0.0003 


-0.0003 


-0.0003 


Sum(//) 


0.0016 


0.0062 


0.0088 



For the range of scales between [i = 0.5 GeV and = 1 GeV the variation is 0.007, or 
27% of the RHS of the DMO sum rule. This variation with scale is of course counterbalanced 
by the p 6 counterterms. However, as just explained, if an estimate for these counterterms is 
used where one has no control over the scale at which the counterterm is fixed, the strong 
variation with scale translates into a large uncertainty in the determination of the LHS 
of the DMO sum rule. Consequently, it will be possible to disentangle L^{M P ) from the 
mass corrections on the LHS of the DMO sum rule only if one can perform an independent 
determination of R^ s — f?y s directly from physical observables. 

To close the discussion, we present some remarks concerning existing and future work in 
this direction. A related process in which L±q contributes at order p 4 is the radiative pion 
decay n — > ev e ^. In this process, however, it is the combination + which occurs, and 
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the structure dependent part of the amplitude is actually used to determine this combination 
of p A counterterms.0 At order p 6 , additional counterterms will contribute, some of them not 
occuring in the two-point functions considered here. Bijnens and Talavera have calculated 
the radiative pion decay amplitude to two- loop order although in chiral SU(2). Using 
the hypotheses of resonance saturation (with just vector- and axialvector resonances), they 
estimate the p 6 counterterm contributions to be small. The sum of the two-loop contributions 
is large, however, leading to an enhancement of about 15% for 2Zg — Zg at the scale M p . The 
advantage of an SU(2) calculation is that mass corrections are always suppressed by factors 
of M%. The large uncertainties due to scale-dependence of resonance saturated counterterms 
found in our SU(3) calculation are thus largely avoided. The comparsion of our results to the 
DMO sum rule in its SU(2) version will thus be an interesting avenue for future work. On the 
other hand, working in chiral SU(3) offers the possibility of studying additional processes 
with external currents/Goldstone Bosons containing strangeness. For instance, there are 
the radiative K i2 decays P2L where more precise data is expected in the near future 



However, because at one-loop order K — > tvfy involves the combination + Ly^ , learning 
more about just will require the study of other processes in which enters separately. 
The meson form factors of vector current matrix elements, (P|V^|Q), where P,Q = tt, K, r], 
offer one possibility to study mass corrections to Lg*\ Partial results on a combination of 



such form factors have been reported recently |24]. It remains to be seen whether the mass 



corrections to the individual form factors can be disentangled from the contributions of the 
low energy constant Lg°\ 



To conclude, we have used our recent two-loop calculation of isospin and hypercharge 
vector and axialvector current propagators in chiral perturbation theory to re-analyze the 
DMO sum rule. The resulting relation, valid to second order in the light quark masses, is 
summarized in Eqs. (|8|)-([T0|). Since we work in SU(3) xSU(3) chiral symmetry, both pion and 
kaon mass dependence is present. A phenomeno logical determination of the DMO spectral 
integral yielded a numerical constraint on the combination of p 4 and p 6 counterterms, given 
in Eq. ([19|) , and prospects for extending the analysis were considered. As regards the overall 
question of how quickly the chiral expansion is converging, we find individual contributions 
at two-loop level to occur at roughly the 30% level. 
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2 In SU(2) notation, the relevant combination of counterterms is 2Z£ — Iq. Its relation to the SU(3) 
couplings is known only to one-loop order. j?]] A direct comparison with our analysis is therefore 
not yet possible. 
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